New insights into the entanglement of disjoint blocks 
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Abstract - We study the entanglement of two disjoint blocks in spin-^ chains obtained by merging 
solvable models, such as XX and quantum Ising models. We focus on the universal quantities that 
can be extracted from the Renyi entropies Sa- The most important information is encoded in 
some functions denoted by Fa. We compute F2 and we show that Fa — 1 and F^.n.^ corresponding 
to the von Neumann entropy, can be negative, in contrast to what observed in all models examined 
so far. An exact relation between the entanglement of disjoint subsystems in the XX model and 
that in a chain embodying two quantum Ising models is a by-product of our investigations. 
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In the last two decades, the entanglement entropy has 
become one of the most important indicators of critical 
behavior in many-body systems [Hll]. In the neighbor- 
hood of a quantum phase transition the area law [3 , i.e. 
the proportionality between the entropy and the contact 
urface with the remainder of the system, does not work in 
ID, where the entropy turns out to be proportional to the 
logarithm of the characteristic subsystem's length. The 
feictor in front of the logarithm is generally universal [3H6] 
and, in conformal systems, is proportional to the central 
charge [7]. In fact, the entanglement entropy provides one 
of the most accurate ways of detecting the value of the 
Central charge. 

. Many recent works [8HT5] have focused on the entangle- 
liient of disjoint intervals in conformal systems, because it 
is also sensitive to universal details of the conformal field 
theory (CFT) different from the central charge. Renyi en- 
tropies = izr^ logTrp" of two disjoint blocks of lengths 
and ^2 at the distance r, as long as a is integer, are 
proportional to the logarithm of the four-point function of 
particular twist fields [10]. For any given a, a (bounded) 
universal function of the four-point ratio x = , ^ , x , 
usually denoted by is the main information, besides 

the central charge c, that can be extracted from the Renyi 
entropy 
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The invariance of the entropies under interchanging 



the subsystem with the rest manifests in the identity 
Pa{x) =Fq,(1 — x). The entanglement of a single spin 
block can be recovered in the limit as r approaches 00, 
and hence x ^ 1. Indeed, when the blocks are far enough, 
the reduced density matrix of the two blocks is expected 
to factorize in the tensor product of the reduced density 
matrices of each block (interaction is local). Thus, the 
Renyi entropies (in particular the entanglement entropy) 
become simply the sum of the entropies of the blocks. This 
means that, normalizing the universal function in such a 
way that Fa(0) = i^a(l) = 1, the constant c'^ in eq. ([T]) is 
exactly the additive constant of the single interval [7] 
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For non- interacting fermions Fo^{x) is equal to 1 for any 
X [16] , however this is peculiar to the fermionic representa- 
tion of non-interacting models [TT] , whilst in the spin rep- 
resentation, as exactly computed in refs. [T0l[l3] in CFT, 
Foc{x) is a complicated function. It has been determined 
only for integral values of a for the free boson compactified 
on a circle [TQ1IT5] and for the quantum Ising model [TSp^ . 
The universal function F^.Ar.(^) = dcxFcx{x)\^_^^ corre- 
sponding to the von Neumann entropy, is still unknown 
for any model. 

One of the goals of this work is to provide exact results 
for Fq,(x), especially for ^2(^)7 in a family of models whose 
continuum limit is conformal. In particular, the finite-size 
scaling can be obtained by substituting any length ^ with 
the chord length ^ sin ^ [7 , where L is the chain's size. 
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Incidentally, in refs. [T7l[T8] it has been shown that the 
excited states associated to primary operators in a CFT 
have instead peculiar finite-size scalings. The results pre- 
sented in this Letter are the first analytic ones obtained 
in a framework different from CFT. 

In light of the numerical investigations [TTl[T3l[T5] and 
of the CFT predictions p!Q[ ri4] , one could wonder whether 
some properties, always observed in Fq,(x), are indeed gen- 
eral. We are referring to the positivity of Fa{x) — 1 and 
of Fy^N^^)^ well as to the concavity of Fy^N^x). The 
latter could be thought as a consequence of some general 
property, in the same way as strong subadditivity implies 
that the von Neumann entropy of a single interval is a con- 
cave function of the block's length [19 ,20 . As a matter of 
fact, some of the models considered here have Fa{x) < 1, 
as well as F^.at. (x) < 0. This is rather surprising since in 
ref. [14] it has been shown that the first order of the short- 
length expansion of Fa{x)^ corresponding to the limit in 
which the distance between the blocks is much larger than 
their lengths, is positive. We will see that the two things 
are in fact compatible. 

The key to our results lies in the non-locality of the 
fermionic mapping of spin-^ chains. As discussed in 
ref. [12], and then explained in detail in ref. [TTj, the en- 
tanglement of disjoint subsystems in the fermionic rep- 
resentation of non-interacting spin-^ chains is different 
from that in the spin representation, because the Jordan- 
Wigner transformation 

4 = n-X' -ith a+ = ^L±M, (3) 

is non-local. Igloi and Peschel in ref. [12] have given a clear 
evidence of this fact by comparing the sub-lattice entangle- 
ment {e.g. the entanglement of the subsystem consisting 
of all even sites) in the two representations: the entan- 
glement entropy, which is extensive in both cases, has a 
different slope. In spite of this, the consequences of the 
inequivalence are often surprising, and the models consid- 
ered here constitute further examples of counterintuitive 
behaviors. 

Finally, studying the entanglement of disjoint subsys- 
tems in a chain embodying two quantum Ising models, we 
find an exact correspondence with the XX model, which 
complements the results obtained by Igloi and Juhasz in 
ref. [21] for the entanglement of spin blocks. 

The model. — We consider a class of spin-^ chains 
closely related to the models introduced by Suzuki in 
ref. [22]. This is an extension of the XY model in which 
the Hamiltonian has a longer-range interaction, and yet it 
reduces to free-fermions after the J-W transformation (|3]). 
The strategy is to "move away" the spin operators adding 
a cr^-string between those that are non-local in the J-W 
fermions (<j^^^^) , so that no string appears in the fermionic 
representation, e.g. 

, ^x ^X f A\ 



We focus on chains that can be constructed by merging 
XX and Ising models, which we call "constitutive models" . 
We consider in particular the Hamiltonians 

XX- : El ^^^Ui • • • <+n-i^r+n + h.c. 

^XX X Is + ^Iz-i^lz4z+i + ^li ' 

where the names (on the left) refer to the constituents: 
the Hamiltonians are sum of XX and Ising models that act 
on different fermionic degrees of freedom; the constitutive 
Hamiltonians are obtained by restricting the sum in eq. ^ 
over the indices that have the same remainder after the 
division by the number n of sub-chains. The models in 
eq. dU can be solved by diagonalizing separately the n 
commuting Hamiltonians. For example, by merging two 
quantum Ising models we eventually get 

^Is2 = XI {^k,e^k,e + bl obk,o - l) , (6) 
k 

where bk,e and bk,o are the Bogolioubov fermions that di- 
agonalize the Ising model, constructed with the even and 
the odd J-W fermions, respectively; Sk = Sk,e = £k,o is the 
dispersion relation. By looking at the low-energy excita- 
tions, we find that the dispersion relation is approximately 
linear and the number of chiral modes {i.e. the number 
of zeros of the dispersion relation above the ground state) 
is equal to the sum of the chiral modes of the constitu- 
tive models. This observation is crucial to infer the value 
of the central charge of the underlying CFT, which turns 
out to be the sum of the central charges of the n con- 
stitutive models. Thus, the long distance physics in the 
examples (|5]) is expected to be described by a CFT with 
central charge n, ^, and |, respectively (c = 1 in the XX 
model and c = ^ in the Ising one). 

We point out that all Hamiltonians with the same struc- 
ture as in eq. ^ are equivalent (up to involution opera- 
tors multiplying the boundary terms) to the corresponding 
Hamiltonians without the cr^-string in eq. (H]), e.g. 

H,,2 ^ Hi, © i/is = X afaf^, + af . (7) 
I 

To avoid confusion we call the latter ones "n-copy Hamil- 
tonians", although in general they could be sum of dif- 
ferent Hamiltonians. The equivalence is realized by the 
transformation 

^nj-\-s^nl-\-s ' (8) 

i<nl-\-s s'<s j j<l 

This can be seen as follows: starting from the n-copy 
Hamiltonian, we apply n J-W transformations (jSj) re- 
stricted to each sub-chain (the last product in eq. ([8])). 
The new Hamiltonian appears like the fermionic represen- 
tation of the merged one, but the operators, which behave 
like fermions inside of their own sub-chain, commute with 
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those of the other sub-chains. We indicate with a^^^ the 
complete -string of the sub-chain s (a^^| commutes with 
H and anticommutes with the fermions in s). In order to 
restore the correct algebra, we multiply the operators in s 
by a^^/| if and only if the operators in 5' are not multiplied 
by a^^i (this is the origin of the central product in eq. (|8|)). 
Finally, the inverse J-W transformation (the first product 
in eq. ([8])) maps the Hamiltonian into the merged one. 

The model obtained by merging n constitutive chains 
and the corresponding n-copy model have different entan- 
glement features. We sum up firstly the properties of the 
n-copy model. 

n-copy model. — Because we are considering spa- 
tial entanglement, the entanglement entropy in the ground 
state of an n-copy Hamiltonian is sensitive to its simple 
structure: the reduced density matrix of any subsystem 
is the tensor product of the reduced density matrices re- 
stricted to the subspaces in which the constitutive Hamil- 
tonians act, e.g. in the model His © ^is? eq. ([7j), we find 
the factorization 



jseis 
Pa 



'p1 
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where (Ao) is the subsystem consisting of the even 
(odd) sites of A, and p^^ is the reduced density matrix 
(RDM) in the Ising model. In particular this means that 
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and 



F„{x) = llFi^\x), 



(11) 



where the index i runs over the constitutive models {of. 
eq. ([1])). The underlying CFT is the tensor product of 
the theories describing the low-energy excitations of the 
constitutive models. No auxiliary information can be ex- 
tracted from the entanglement in the ground state of n- 
copy Hamiltonians, apart from those already known from 
the constitutive models. 

Two disjoint blocks. — We now come back to our 
model. Because the entanglement of a single interval is 
equal both for spins and J-W fermions (the J-W transfor- 
mation, although non local, mixes only degrees of freedom 
inside of the block), eq. (p!Q|) continues to be valid: we do 
not find new information studying the block entanglement. 
However, in chains with periodic boundary conditions, the 
RDM of disjoint blocks is different from the corresponding 
density matrix in the fermionic representation. In ref. |TT 
it has been shown that the RDM of two disjoint blocks 
A = Ai U A2 (we call Bi and B2 the blocks between Ai 
and A2) is equivalent to the sum of four operators (two 
genuine density matrices and two "fake" ones) 
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Pa +<2PA «A2 



z (Bi) z 



(12) 



which for non-interacting models are exponentials of 
quadratic forms. The first operator is the fermionic den- 
sity matrix. We indicated with the a^-string, product 
of every in R. The fake density matrix is defined as 



3.)P 



(Bi) 



IVB,uBj|*o){*o|a|J 



(13) 



where |^o) is the ground state. In fact, the main dif- 
ference from the n-copy model is that the factorization, 
shown in the example ([9]), holds now for each operator in 
eq. (p!2|) separately (in particular the expectation value of 
the string factorizes). The second moment TTp\ in a chain 
constructed with n identical constitutive models is given 
by (see ref. [11 for the basic case n = 1) 



{r?}" + {ri,r2r 



2„ {ri}"-{r3,r4}' 



_ (14) 

where we used the same notation as in ref. [11^ for the 
correlation matrices: Fi = I — (a a) is the fermionic 
correlation matrix, with the Majorana fermions defined 
as a2i-i = c] + Ci and a2i = i{ci - c-); F2 = P2F1P2 and 
F4 = P2TsP2^ where P2 inverts the sign of the components 
corresponding to the second block A2; F3 is given by 



r3 = Fi 



^ABi^BiBi^BiA 



(15) 



where the double subscripts take into account restrictions 
to rectangular correlation matrices, i.e. the first (second) 
subscript identifies the region where the row (column) in- 
dex runs {e.g. Taa = Fi). {F,F'} is the product of 
the eigenvalues of (I + FF')/2 with halved degeneracy (see 
ref. fll, for further details). We stress that in the correla- 
tion matrices above each length is divided by n, indeed the 
fermionic space of the constitutive models is the n^^ root 
of the total space. The generalization to distinct constitu- 
tive models is straightforward: we must substitute powers 
with products in eq. ([14]). 

In order to get the universal function ^2(^)7 because 
{F^} = {Fi,Fi} is the fermionic analogous of eq. ([H]), 
corresponding to F2{x) = 1 (free fermions), we can divide 
Trp^ by {Tf}^^ which takes into account the actual value 
of the central charge c and of the constant C2 {cf. eq. ([1]) 
and eq. 1^) 



^(n), , 1 , {ri,F2r + (a|^)^-({Fir-{F3,F4n 
^2 (^) = 2+ 



2\n 



2{r?} 

(16) 

It is remarkable that eq. ([16]) gives direct access to the uni- 
versal function, which is a subleading term of the Renyi 
entropy 5*2 {cf. eq. ([1])). We point out that the knowledge 
of F2^\x) for n = 1,2,3 is sufficient to get the function 
for any n. However, if we are interested in Fq^^\x) (or in 
F^^\x)) we need 4 (or 18) values of n to compute the uni- 
versal function for any n, because of the increasing num- 
ber of terms in the formulae analogous to eq. ([T^ (see eqs. 
(58,59) of ref. [II]). 
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Fig. 1: The terms in eq. {[T9|) for ^ + r = 25, 50, 100; symbols 
are sorted by increasing opacity and decreasing size. In the 
legend a = {Fi, r2}/{rf }, b = {a%^f {Ti}/{Ti}, and c = 
(o^Si)^ {^3, r4}/{r?}. Data converge quickly to the functions 
proposed in eq. (|19p (solid lines). 



Equations (p!2]) . (p!4|) , and ([T6]) are general: they hold 
for any set of non-interacting constitutive models. Equa- 
tion (p!2]) is valid also in the presence of interaction. We 
specialize now the formulae to XX and Ising sub-chains. 

XX model. — The entanglement entropy of two dis- 
joint blocks in the XX model has been investigated in rel- 
atively small chains by Furukawa et al. in ref. [9j by exact 
diagonalization techniques. In ref. [11] the CFT predic- 
tion IHITO] 



{x) 



1 + yr 



(17) 



has been checked numerically in the thermodynamic limit. 
Now, in order to find the analogous results for n > 1, i.e. 
for the ground state of the Hamiltonian 



l+n-l^l^rn 



h.c. , 



(18) 



we should determine the behavior of each of the terms 
in eq. ([H]). The correlation matrix Fi is a Toeplitz ma- 
trix in which a block of rows (associated to the region 
Bi) and the corresponding block of columns have been 
removed. Therefore, the basic problem of computing the 
Renyi entropies of disjoint blocks of the same length in the 
fermionic representation of XX chains can be traced back 
to the asymptotic behavior of the determinant of large 
block- Toeplitz matrices with Fisher-Hartwig singularities. 
However, we are not aware of useful mathematical theo- 
rems or conjectured for this kind of matrices. In addition, 
the other matrices appearing in eq. ([H]) have structures 
even more complicated. 

In any case, we provide numerical evidence that each of 
the terms summed in eq. ([16]) is a simple function of x, 
namely 



{ri,r2} 



(19) 
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Fig. 2: The universal function F2{x) for the XX"^ model 
(eq. ([T8)) with n = 3). The universal function is evidently 
smaller than 1. The continuous curve is eq. (1201) with n = 6. 



while the last term of eq. ([T6|) is subleading. In fig. [T] the 
proposed scaling functions are checked against numerics in 
the thermodynamic limit for various subsystem's configu- 
rations with blocks of the same length ii = £2 = £• The 
agreement is excellent. In general, by inserting eq. ([19]) 
into eq. ([T6]) we get 



1 + (1 - a;)t +; 



(20) 



For n = 2 the function simplifies to F2 {x) = 1. However, 
for n > 2 the universal function becomes smaller than 1, 
and in the limit n ^ 00 we get F2^^(x) = | for any 
X G (0, 1). In fig. [2] we report the universal function cor- 
responding to n = 3 for two values of ^ + r. 

Equations (jl6|2Qp are the main results of this Letter. 
In particular eq. ([20]) is the first expression for the univer- 
sal function F2 (x) obtained in a framework different from 
CFT. Furthermore, for n > 2 it is the first example in 
which F2{x) < 1. 

Calabrese et al. in ref. [14 have determined the short- 
length expansion for the free boson and Ising model 



F2{x) = l 



(21) 



where y is the scaling dimension of an operator coming 
from the operator product expansion of twist fields, and A/" 
is the number of inequivalent correlation functions given 
the same contribution. By series expanding eq. ([20]) for 
small X we find 



1 



(22) 



As observed in ref. [14], when there is no operator in the 
theory whose correlation function can contribute to 0(x), 
the linear term in eq. ([2T]) must cancel with the 0{x) term 
coming from the function multiplying Fa{x) in Trp", i.e. 
it must be equal to — This is exactly what hap- 

pens here, being the central charge c = n {of. eq. ([TO])). 
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If eq. (|2T]) can be applied, by direct comparison between 
eq. (|22|) and eq. (|2l]) we recognize y = ^ and Af = 4'^/2. 

The fact that F2^\x) < 1 for n > 2 is essentially due to 



the 0{x) contribution in eq. 

Incidentally, the short-length expansion in the n-copy 
model is characterized by the same exponent y = ^ of the 
XX model and JV = 2n, being the n-copy model sum of n 
identical models. 

From the analytic continuation of the short-length ex- 
pansion of Fq,(x), it can be deduced the expansion of 



F..iv.(^)=A/'(f) 



x\y^/^T{y^l) 



4r(^ + f) 3 

By substituting the value of y we get 

„v^r(^) c 



F(l(x) = (4a:) 



8r(^) 



(23) 



(24) 



In particular the universal function F^.Ar.(^) vanishes for 
n = 2. For n > 2 F^^^ (x) is negative (at least for 
small enough x), in contrast to what observed up to this 
time |9j l3, 23]. We point out that if i^v.Ar.(^) is negative 
then the first order of the small-length expansion must be 
linear. This is compatible with the constraints given by 
strong subadditivity [20] . 

The models obtained by merging XX chains, eq. (fTS)) . 
have a further property: they commute for different val- 
ues of n, and in particular they commute with the XX 
Hamiltonian (n = 1). Thus, from eq. (|2Ql) and eq. (plQ)) 
with C(^) = 1 and {02) (i) = C2 the additive constant for 
the block Renyi entropy ^2 (see ref. [24] for the analytic 
expression), we get the entropy of two disjoint intervals in 
some excited states of the XX model. In this way we have 
extended the analysis of refs. [23l[25l[26], focused on the 
entanglement of spin blocks. 



Ising model. 

Ising chains 



We now consider the merging of n 



His 



2^ ^i+i 



I 



(25) 



The first approximate numerical results for the entangle- 
ment entropy of two disjoint blocks in the Ising model 
(n = 1) have been obtained for a = 2 in ref. |i3j by using 
a tree tensor network algorithm and Monte Carlo simula- 
tions. In ref. |11J the CFT prediction 



(26) 



has been checked against exact numerical data. As done 
before for the XX model, in order to get the universal 
function ^2"^^ (x) we must determine the behavior of each 
of the terms in eq. (p!6|) . The result, as shown in fig.[3l has 
the same structure as in the XX model {of. eq. (plO]) ) 




Fig. 3: The terms in eq. (O for ^ + r = 25, 50, 100; symbols 
are sorted by increasing opacity and decreasing size, a, 6, and c 
are the same quantities defined in fig. [1] Data converge quickly 
to the functions proposed in eq. (|27)) (solid lines). 



and the last term of eq. ([T6]) is again subleading. The uni- 
versal function corresponding to merging n Ising models 
is then given by 



1 + (1— X)4 + 



(28) 



This formula is identical to eq. (|2Q|) . provided that n is 
substituted by ^. In particular the central charge c and 
the exponent y are both equal to ^. For n = 2, namely 



'l ^l-\-l^l-\-2 



(29) 



both the additive constants of the single interval {of. 
eq. dHJ) and the universal function F2{x) are equal to those 
corresponding to the XX model. In fig. |l]we report both 
F2 {x) and F3 (x), showing agreement with the CFT 
predictions relative to the XX model. We expect this re- 
mains true for any a {of. eq. ([T])). The relation between 
the entanglement in the XY model and that in the Ising 
model has been already discussed in ref. [21] , where the au- 
thors focused on the entanglement entropies of spin blocks. 
Here we have established a correspondence between the 
entanglement of disjoint blocks in the XX model and that 
in two merged Ising chains, eq. (|29|) . The two models 
display different corrections to the scaling, however the 
asymptotic behavior of Renyi entropies coincides. 

In contrast to the XX model, the Hamiltonians ([25]) do 
not commute with the Ising Hamiltonian. 

The third example in eq. (|5]) shows the merging of an XX 
chain with an Ising one. The Hamiltonian is not transla- 
tional invariant, however the symmetry will be eventually 
recovered in the continuum limit, as it happens for the 2- 
copy Hamiltonian. Equation ([16]) still holds for this kind 
of models and in general we find 
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Fig. 4: The universal functions F^^^ (x) and F^^^ (x) ioi l + r = 
25, 50, 100 in the Ising^ chain (eq. (|25)) with n = 2); symbols are 
sorted by increasing opacity and decreasing size. Data converge 
to the universal functions of the XX model (continuous curves). 



where nxx and njs are the numbers of XX and Ising consti- 
tutive models, respectively. Because of the correspondence 
between XX and merged Ising chains, this result does not 
add significant information. However, in general, merging 
different chains results in universal functions that are not 
trivially related to the constitutive ones. 

Summary and discussion. — We have investigated 
the entanglement entropy of two disjoint blocks in exactly 
solvable spin-^ chains. We derived the universal function 
F2{x) for the XX and the Ising model, reproducing the 
known results from CFT. In addition, we computed F2{x) 
for an entire class of chains that embody models with a 
free-fermion representation. We have shown that Fo^{x) 
can be smaller than 1, as well as Fy^N.ix) negative, by 
providing examples in which this indeed happens. Finally, 
we revealed a correspondence between the XX model and 
two merged Ising chains, extending the analysis of ref. [21] 
to the entanglement of disjoint subsystems. 

In this paper we considered a = 2, however the analysis 
of Renyi entropies with a > 2 could be useful to under- 
stand the connection between the CFT formalism and the 
method of ref. pT], based on free-fermion techniques, per- 
haps making progress toward the problem of computing 
the universal function F^.at. (^) of the von Neumann en- 
tropy. 

In light of the existence of models in which F^. at. {x) < 0, 
it becomes worth determining the lower bound of the 
universal function, for example investigating the conse- 
quences of strong subadditivity [20] . 

Finally, we have seen that the XX Hamiltonian com- 
mutes with the models embodying XX chains. However, it 
does not commute with the corresponding n-copy Hamil- 
tonians. Thus, it could be interesting to study the features 
of the stationary state that describes the system late times 
after the quench resulting from joining together the n de- 
coupled XX sub-chains to form an XX chain. 
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